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Amplification of deterministic disturbances in inertialess shear-driven channel flows of 
viscoelastic fluids is examined by analyzing the frequency responses from spatio-temporal 
body forces to the velocity and polymer stress fluctuations. In strongly elastic flows, we 
show that disturbances with large streamwise length scales may be significantly amplified 
even in the absence of inertia. For fluctuations without streamwise variations, we derive 
explicit analytical expressions for the dependence of the worst-case amplification (from 
different forcing to different velocity and polymer stress components) on the Weissenberg 
number (We), the maximum extensibility of the polymer chains (L), the viscosity ratio, 
and the spanwise wavenumber. For the Oldroyd-B model, the amplification of the most 
energetic components of velocity and polymer stress fields scales as We^ and We^. On 
the other hand, finite extensibility of polymer molecules limits the largest achievable 
amplification even in flows with infinitely large Weissenberg numbers: in the presence of 
wall-normal and spanwise forces the amplification of the streamwise velocity and polymer 
stress fluctuations is bounded by quadratic and quart ic functions of L. This high am- 
plification signals low robustness to modeling imperfections of inertialess channel flows 
of viscoelastic fluids. The underlying physical mechanism involves interactions of poly- 
mer stress fluctuations with a base shear, and it represents a close analog of the lift-up 
mechanism that initiates a bypass transition in inert ial flows of Newtonian fluids. 



1. Introduction 

Newtonian fluids transition to turbulence under the influence of inertia. In stark con- 
trast, recent experiments have shown that flows of viscoelastic fluids may undergo a 
transition to a time-dependent disordered flow state and become turbulent even when 



inertial forces are considerably weaker than viscous forces (Larson 2000 Groisman 



Steinberg 2000 2004 Arratia, Thomas, Diorio & Gollub 2006). Since viscoelastic fluid 



flows are often encountered in commercially important settings, understanding transition 
to elastic turbulence in such flows is important from both fundamental and technologi- 
cal standpoints. In polymer processing, for example, elastic turbulence is not desirable 



because it compromises quality of the final product (Larson 1992). But in microfluidic 



devices elastic turbulence can help promote transport, thereby improving the quality of 
mixms dGroisman fc Steinberg||2001[ [Ottino fc Wiggins||2004D. 



Transition in the experiments of Groisman & Steinberg ( 2000 , 2004 ) is thought to be 



initiated by the occurrence of a linear instability that arises from the presence of curved 
streamlines (Larson, Shaqfeh & Muller 1990; Larson 1992). However, the question of 



whether and how transition can occur in channel flows of viscoelastic fluids with straight 
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streamlines remains wide open. Standard modal stability analysis of the upper convected 
Maxwell and Oldroyd-B constitutive equations shows that these flows are linearly stable 
when inertial effects are negligible; yet, they exhibit complex dynamical responses in 
strongly elastic regimes (Yesilata'2002; '2009! [Bonn et a/.||2011[ [Pan et a/.|[2Qll] ). Thus, 
if an inertialess transition can indeed be described using such basic constitutive mod- 
els, it would likely involve finite- amplitude disturbances that would trigger nonlinear 
effects ( [Meu lenbroek et al. 2004 Morozov & van Saarloos 2005). However, the lack of 
a modal instability does not rule out the possibility that the early stages of transition 
can be described by the linearized equations. If non-modal growth is present, initially 
small-amplitude disturbances could grow to a finite amplitude at intermediate times be- 
fore decaying at long times. For sufficiently large disturbance amplitudes the flow could 
enter a regime where nonlinear interactions are no longer negligible. This can induce 
secondary amplification and instability of the flow structures that are selected by the 
linearized dynamics and promote eventual transition to elastic turbulence. 



Hoda, Jovanovic & Kumar (2008, 2009) recently employed tools from linear systems 



theory to study the amplification of stochastic spatio-temporal body forces in plane 
Couette and Poiseuille flows of viscoelastic fluids with nonzero Reynolds numbers. In 



strongly elastic flows, the results of Hoda et al. (2008 2009) indicate that significant 



amplification of streamwise-constant velocity fluctuations can occur even when inertial 
forces are weak. As in Newtonian fluids, this amplification is caused by non-normality of 
the underlying operators and it cannot be predicted via standard linear stability analy- 
sis. Furthermore, recent work of j Jovanovic & Kumar (2010 2011) shows that this large 
amplification arises from the interactions between the polymer stress fluctuations in the 
wall-normal/spanwise plane with the base shear. Through these interactions weak stream- 
wise vortices induce a viscoelastic analogue of the lift-up mechanism which is responsible 
for the creation of alternating regions of high and low streamwise velocities (relative to 
the mean flow), po^vanovic fc Kumar" (2010, 2011) demonstrated significant conceptual 
similarities between this purely elastic mechanism and the well-known inertial vortex 
tilting mechanism that initiates a bypass transition in shear flows of Newtonian fluids. 

Despite this recent progress, analytical results that quantify influence of finite exten- 
sibility of polymer molecules on the amplification of disturbances in channel flows of 
viscoelastic fluids without inertia are still lacking. Such results may provide physical in- 
sight into the early stages of transition and help benchmark direct numerical simulations. 
Analogous results have been extremely helpful in understanding the early stages of tran- 
sition to turbulence in wall-bounded shear flows of Newtonian fluids (Farrell 



loannou 



1993| [Trefethen erZ]|1993| [Jovanovic fc Bamieh||2005[ |Schmid|[2QQ7l ). 



It is worth noting that the problem of determining the amplification of white-in-time 
stochastic forcing is ill-posed when inertia is completely absent. This restricts the results 



of Hoda et al. ( 2008 2009 ) to cases where the flow has finite inertia. Jovanovic fc Kumar 



(2011) used singular perturbation methods to identify the spatial structure of velocity 
and polymer stress fluctuations that exhibit the highest amplification in stochastically 
forced weakly inertial channel flows of viscoelastic fluids. As the influence of inertial 
forces vanishes, it was demonstrated that the velocity fluctuations become white-in-time. 



thereby exhibiting infinite variance. Furthermore, Hoda et al. (2008, 2009) and Jovanovic 



fc Kumar (2010, 2011) employed the Oldroyd-B constitutive model to investigate energy 



amplification of velocity and polymer stress fluctuations. However, since the Oldroyd-B 
model allows the polymers to stretch indefinitely, examining the role of finite extensibility 
of polymer molecules on the amplification of disturbances remains an open question. 

In the present work, we address these issues by examining the worst-case amplifica- 
tion of deterministic disturbances in inertialess (i.e., creeping) shear-driven channel flows 
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of viscoelastic fluids. We consider spatially distributed and temporally varying forcing 
that is purely harmonic in the horizontal directions and time, and deterministic in the 
wall-normal direction. The motivation for studying creeping flows arises from the obser- 
vation that viscoelastic fluids can become turbulent even in low inertial regimes, i.e., at 
small Reynolds number (Larson 1992 Groisman fc Steinberg^ ^2000 J . Furthermore, the 
present analysis uses the finitely extensible nonlinear elastic Chilcott-Rallison (FENE- 
CR) model (Chilcott & Rallison 1988"), which captures the finite extensibility of the 
polymer molecules. It is well-known that, for infinitely extensible polymers, the FENE- 
CR model simplifies to the Oldroyd-B model. 



With our approach, we show that velocity and polymer stress fluctuations can ex- 
hibit significant amplification and that the most energetic flow structures have large 
streamwise length scales. In the absence of streamwise variations, we derive explicit ex- 
pressions for the worst-case amplification from different components of the forcing field 
to different components of velocity and polymer stress fluctuations. For the Oldroyd-B 
model, the wall-normal and spanwise forces induce amplification of the streamwise com- 
ponents of velocity and polymer stress fields that scales quadratically and quarticly with 
the Weissenberg number. On the other hand, we demonstrate that finite extensibility of 
the polymer molecules saturates the largest achievable amplification even for flows with 
infinitely large Weissenberg numbers. The functions that bound the worst-case amplifi- 
cation of the velocity and polymer stress fluctuations scale quadratically and quarticly 
with the largest extensibility of the polymer molecules. We also identify the viscoelastic 
analogue of the well-known inertial lift-up mechanism as the primary driving force for 
high flow sensitivity; the underlying mechanism arises from interactions of polymer stress 
fluctuations with a base shear and it is facilitated by spanwise variations in flow fluc- 

). This non-modal amplification may provide a route 
by which infinitesimal disturbances can grow to finite amplitude and perhaps trigger a 
transition to elastic turbulence in channel flows of viscoelastic fluids. 



tuations (Jovanovic fc Kumar||2Qll 



To facilitate development of explicit analytical expressions for worst-case amplification 
of velocity and polymer stress fluctuations, we restrict our study to an inertialess shear- 
driven channel (Couette) flow of FENE-CR fluids. Even though the current analysis can 
be readily applied to the FENE-P model and to a pressure-driven channel (Poiseuille) 
flow, more complicated base state removes algebraic convenience encountered in Couette 
flow of FENE-CR fluids. We note that all physical mechanisms identified in this paper 
remain at play in inertialess Poiseuille flow of both FENE-CR and FENE-P fluids. 



The rest of this paper is organized as follows: in § [2j we use a simple example to 
illustrate how techniques from control theory can be used to quantify amplification of 
disturbances and robustness to modeling imperfections. In §[3j we describe the governing 
equations for inertialess channel flows of FENE-CR fluids, provide the evolution model, 
and briefly discuss the essential features of the frequency response analysis. In § |4j we 
examine the frequency responses of three-dimensional velocity fluctuations in inertialess 
Couette of FENE-CR fluids. In § [Sj we provide analytical expressions for the worst- 
case amplification from the forcing to velocity fluctuations using a streamwise-constant 
linearized model. We also identify the spatial structures of the dominant forcing and 
velocity components and demonstrate the importance of the viscoelastic lift-up mecha- 
nism. In § |6) we study the dynamics of streamwise-constant polymer stress fluctuations. 
Finally, in § [7| we summarize the major contributions of this work and discuss future 
research directions. 
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Figure 1. Block diagrams of (a) system ( |2.1| ); and (b) system (2.1) connected in feedback with 
norm-bounded unstructured uncertainty 1. Here, s G C denotes the temporal Laplace transform 
variable, and the transfer function, from d to cp, is determined by H{s) = R/ ((s + Ai)(s + A2)). 



2. The role of uncertainty: an illustrative example 

In the course of addressing the issue of disturbance amplification, the present paper 
provides insight into the robustness of viscoelastic flow models. In this section, we briefly 
summarize how the tools from control theory facilitate quantification of sensitivity and 
robustness of a system to modeling imperfections via frequency response analysis. The 
approach taken in the present paper is closely related to classical frequency response 
studies of systems arising in physics and engineering. For example, in the design of oper- 
ational amplifiers it is well-known that caution must be exercised with models that show 
a low degree of robustness because small modeling errors could cause otherwise stable dy- 
namics to become unstable. Similar ideas have found use in fluid mechanics, including the 
analysis of the early stages of transition in shear flows of Newtonian fluids ( Schmid||2QQ7 ). 



To fix ideas, let us begin with a simple system of two coupled first-order differential 
equations 



" 01 ■ 




' -Ai 







" 4>i 


4>2 
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-A2 _ 







f = [0 1 ] 



(2.1) 



where and 02 are the states, d is the disturbance, and (f is the output. We assume 



positivity of scalars Ai and A2, which guarantees modal stability of (2.1). Equivalently, 
this system can be represented via its block-diagram in figure |l(a) Clearly, we have 



a cascade connection of two stable first-order systems with parameter R determining 
the strength of coupling between the two subsystems. For Ai ^ A2 the solution of the 
unforced problem, i.e. with (i = 0, is determined by 



hit) 
hit) 



-Alt 



2(0) 



R 



A2 — Ai 



(2.2) 



Thus, the initial conditions in (pi create monotonically decaying responses of with a 
rate of decay determined by A^. In contrast, the response of (/)2 arising from the initial 
condition in (pi is characterized by two competing exponentials, and it vanishes both for 
t = and for asymptotically large times. For finite times, however, transient growth, 
directly proportional to the coupling coefficient i?, is exhibited with the largest value of 
transient response taking place at t = (l/(Ai — A2)) log (A1/A2). This transient growth 
does not require the presence of near resonances (i.e., Ai ~ A2) or modes with alge- 
braic growth (i.e., Ai = A2); it is instead caused by the non- normality of the dynamical 
generator in \2.1\ . 

We note that transient growth represents one particular manifestation of the non- 
normality of the dynamical generator in the above example (see Grossmann ( 2000 ) for a 
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(can be nonlinear or time-varying) 

Figure 2. Block diagram of a system connected in feedback with norm-bounded unstructured 

uncertainty F. 



comprehensive treatment). Additional features can be observed by analyzing (2.1) in the 
frequency domain. The frequency response is obtained by evaluating the transfer function 
H{s) (from input d to output cp^ (p{s) = H{s)d{s)) on the iuj axis, where G M is the 
temporal frequency and i = ^/^. The largest value of [^^(ico')! determines the so-called 
Hoo norm. This measure of input-output amplification has an appealing physical inter- 
pretation: it quantifies the worst-case amplification of finite-energy disturbances ( |Zhou| 
et a/.|p!996 ). In the above example, ||i^||oo = |^|/(^i^2) indicates the existence of a unit 



energy disturbance that generates output whose energy is given by |i^P/(AiA2)^. 

Furthermore, the Hoo norm has an interesting robustness interpretation which is closely 
related to the analysis of pseudospectra of linear operators ([ Trefethen fc Embree|[2QQ5) . 
Namely, ||i^||oo determines the size of modeling uncertainty, d{s) = r{s)ip{s), that can 
destabilize the system; see figure [2] for an illustration. This uncertainty may arise from 
the inevitable imperfections in the laboratory environment or from the approximate na- 
ture of the governing equations (caused by, e.g., high-frequency unmodeled dynamics, 
parametric variations, neglected nonlinear it ies, or crude physical assumptions made in 
modeling). In particular, system ( |2.1|) with d{s) = T{s)(f{s) can be represented by a 
feedback interconnection in figure l(b)[ If, apart from being norm-bounded, there are 
no structural restrictions on uncertainty F, then the necessary and sufficient condi- 
tion for stability of a feedback interconnection in figure |l(b)| is given by the so-called 
small-gain theorem, ||r||oo < l/||i^||oo. In the above example, this condition simplifies 
to ||r||oo < AiA2/|i^|. In particular, it is easy to establish the existence of a constant 
gain uncertainty, F(s) = 7 = const., of magnitude larger than AiA2/|i^| that makes the 
system in figure 1(b) [ unstable. 

The above example illustrates that, in systems with non- normal generators, the eigen- 
values may represent misleading measures of both transient growth and input-output 
amplification. While they successfully predict system's behavior for asymptotically large 
times, they may fail to capture important aspects of short-time behavior, disturbance 
propagation, and robustness. In particular, the coupling R between subsystems in (2.1) 
plays a crucial role in determining transient and input-output features of the system's 
response: large values of R signal large transient responses, poor stability margins, and 
large amplification of disturbances. In the absence of the coupling, i.e., for = 0, the 
stability margins of subsystems in (2.1) are determined by Ai and A2. In contrast, for 
non-zero R this margin is determined by AiA2/|i^|, thereby indicating that, even for sub- 
systems with large stability margins, small uncertainties can have a destabilizing effect 
on the overall system if the coupling between the subsystems is large enough. 

The main ideas from the above example extend to multivariable and infinite dimen- 
sional systems. For these problems, the singular values of the frequency response operator 
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can be used to determine input-output amplification in the presence of disturbances. Fur- 
thermore, the analysis of spatio-temporal frequency responses for spatially distributed 
systems can be used to identify prevalent spatial length scales and spatio-temporal pat- 
terns that are most amplified by the syste m's d ynamics. For example, if instead of being 
constant scalars, parameters and R in (2.1) are given by {A^ = + 6^/^^, R = ci/^}, 
then (2.1) can be interpreted as an equivalent of the following system 



(t>2t{x,t) = C(l)ix{x,t) + h2(t)2xx{x,t) - a2(t)2{x,t), 

(p{x,t) = (j)2{x,t), X e 

in the spatial frequency domain after applying the spatial Fourier transform to the 
above system. Here, z^: G M denotes the spatial wavenumber, and (a^,6^,c) denote pos- 
itive reaction, diffusion, and convection coefficients. The /^-parameterized Hoo norm, 
c\K\/{{ai + hin^){a2 + ^2^^)), disappears for = and as z^: ^ 00, thereby achiev- 
ing its peak for non-zero R. This value of n identifies the spatial length scale, 27r/^, 
that has the smallest stability margin and that is most amplified by deterministic distur- 
bances. Thus, convective coupling in reaction-diffusion systems can provide dynamical 
responses that cannot be inferred by analyzing subsystems in isolation. 

We finally note that this simple example captures the essential features of nonmodal 
amplification in wall-bounded shear fiows of both Newtonian and viscoelastic fiuids. In 
Newtonian fiuids the subsystems in figure [T(a)] would correspond to the Orr-Sommerfeld 
and Squire equations and the coupling between them would represent the vortex tilt- 



ing term whose strength is directly proportional to the Reynolds number (Jovanovic & 
Bamieh 2005). In a study focusing on transient growth in inertialess channel fiows of 
viscoelastic fiuids, Jovanovic & Kumar (2010) showed that polymer stretching and the 



Weissenberg number effectively take the role that vortex tilting and the Reynolds number 
play in inertial flows of Newtonian fluids. 



3. Problem formulation 

In this section, we present the governing equations for inertialess shear-driven channel 
flow of viscoelastic fluids. We show how the linearized equations can be cast into an 
evolution form that is amenable to both analytical and computational developments. We 
then provide a brief description of frequency responses and input-output norms, along 
with the numerical tools for computing them. 

3.1. Governing equations 

The non-dimensional momentum, continuity, and constitutive equations for an incom- 
pressible shear-driven channel flow of viscoelastic fluids, with geometry shown in flgurejsj 
are given by ( [Bird et a/.|[l987[ |Larson|[T999| 

/3V^V - VP - i^eV- VV), 



Re\ 

R 



We {{I 
V - V, 



/3) V-T 



We 



VV + (R- VV)-' - V- VR 



(3.1a) 
(3.16) 

(3.1c) 



Here, the overdot denotes a partial derivative with respect to time t, V is the velocity 
vector, P is the pressure, T is the polymer stress te nsor, R is the conformation tensor, 
V is the gradient, and is the Laplacian. System (3.1) governs the behavior of dilute 
polymer solutions with fluid density p, and it has been obtained by scaling length with 
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the channel half- height /i, time with the fluid relaxation time A, velocity with the largest 
base flow velocity Uq, polymer stresses with r]pUo/h^ and pressure with {r]s + Vp) Uo/h, 
where r]s and r]p are the solvent and polymer viscosities. This scaling leads to three pa- 
rameters that characterize the properties of (3.1): the viscosity ratio, /3 = r]s/ {r]s + r]p)] 
the Weissenberg number, We = XUo/h; and the Reynolds number. Re = pUoh/{r]s + 
While the Reynolds number quantifies the ratio of inert ial to viscous forces, the Weis- 
senberg number determines the ratio of the fluid relaxation time A to the characteristic 
flow time H/Uq. 



The momentum (3.1a) and continuity (3.16) equations describe the motion of an in- 



compressible viscoelastic fluid. For given T, the pressure adjusts itself so that the velocity 
satisfies the continuity equation (3.16). In our previous work ( Hoda et a/.|2QQ8 2009 Jo- 
vanovic fc Kumar||2010 , 2011), we used the Oldroyd-B model, which is based on a linear 
bead-spring dumbbell, to relate the polymeric stress tensor to the conformation tensor. 
However, it is well-known that the Oldroyd-B model does not account for the finite exten- 
sibility of the polymer chains. In this work, we address this issue by using the FENE-CR 
model, which utilizes a nonlinear relationship between the polymeric stress tensor T and 
the conformation tensor R ( Bird et aL||1987 ), 



(3.2) 



Here, I is the unit tensor, and the function / (that quantifies the influence of the nonlinear 
spring) is determined by the trace of the conformation tensor, trace (R), and the square 
of the maximum extensibility of polymer chains. 



/ 



L2 



i2 - trace (R) ' 



(3.3) 



Note that R and are made dimensionless with respect to /cT/c, where /c, T, and 
c denote the Boltzmann constant, the absolute temperature, and the spring constant 
of the Hookean dumbbell, respectively. In the limiting case L ^ oo, we have / ^ 1; 
consequently, the nonlinear spring becomes linear and the FENE-CR model simplifies to 
the Oldroyd-B model. 



In a shear-driven channel flow, system of equations (3.1) - (3.2) exhibits the following 
steady-state solution 



R 



[ U{y) 

Rii R21 Ri3 

R12 R22 R23 

Ri3 R23 R33 



l + 2WeV/2 We If 
We/f 1 

1 
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where 



U{y) 



y, 




2 I * ' V * ' i2 

We also note that the first normal stress difference in Couette flow is determined by 

Ni = Ru - R22 = 2{We/f)\ 

In §[5] and §|6]we will show that this parameter, that can take values between and Z^, 
plays a key role in the dynamics of velocity and polymer stress fluctuations. 

In the absence of inertia, i.e. in flows with Re = 0, the dynamics of infinitesimal 
velocity, v, pressure, p, and conformation tensor, r, fluctuations around the base flow 
(v, R) are governed by 



- Vp + (1 - /3) V-T + (3V\ 
V-v, 




(3.4a) 
(3.46) 

(3.4c) 
{3Ad) 



Here, (3 Ad) establishes a relation between polymer stress and conformation tensor fluc- 



tuations. Furthermore, and w are the components of the velocity fluctuation vector 
\- = [ u V w ] in the streamwise (x), wall- normal (?/), and spanwise (z) directions, 
respectively. The momentum equation (3.4a) is driven by a spatially distributed and tem- 
porally varying body forcing, d = [ di ^2 where di, (^2, and ds are the forcing 
fluctuations in the streamwise, wall-normal, and spanwise directions, respectively. In prior 
work using the Oldroyd-B constitutive equations (Hoda et a/.| |2QQ8 



2009 Jovanovic 



Kumar poll ), the three-dimensional body forcing varies harmonically in the horizontal 



directions and stochastically in the wall-normal direction and in time. However, given 



the static-in-time momentum equation (3.4a), white-in-time stochastic disturbances d 



induce white-in-time velocity fluctuations v, and the problem of variance amplification 
(in the absence of inertia) becomes ill-posed ( Jovanovic fc Kumar||20lT ). Hence, in this 
work, we consider the body forcing d to be purely harmonic in the horizontal directions 
and time, and deterministic in the wall-normal direction, and we study the worst-case 
amplification of deterministic disturbances. 



3.2. Model in the evolution form 

We note that (3.4a) and ( |3.46| ) can be simplified by expressing the velocity fields in terms 
of the wall-normal velocity (v) and vorticity {r] = dzU — dxw) fluctuations. This is done 



by first taking the divergence of (3.4a) and using (3.46) to get an expression for p. The 



equation for v is then obtained by eliminating p from (3.4a). The equation for rf can 



be obtained by taking the curl of (3.4a). Finally, by rearranging the components of the 
conformation tensor into the vector 



Worst- case amplification in inertialess Couette flow of viscoelastic fluids 



with 

01 = [ ^22 ^23 ^33]^, 02 = [ ^13 

and by applying the Fourier transform in the x- and 2:-directions, we arrive at the fol- 
lowing static-in-time expressions for v and rj in terms of the conformation tensor and 
body-forcing fluctuations 

V = + d, 
T] = + d. 
In addition, equation (3.4c) can be brought to the following form 

01 = Fii 01 + Fi^ V + Fi^T^, 

02 = F21 01 + F22 02 + F2^ V + F2r7 V- 



(3.5) 



(3.6) 



The operators in (3.5) and (3.6) are defined in Appendix [A| For notational conve- 
nience, we have suppressed the dependence of {v,77,0^,d} on (k, t; /3, We, L), where 
K = {kx, kz) with kx and kz denoting the horizontal wavenumbers. 

The boundary conditions on the wall-normal velocity and vorticity are dictated by the 
no-slip and no-penetration requirements 

v{f^^y = ±1, t) dyV (k, y = ±l,t) = t] (k, y ±1, t) 0. 

We note that there are no boundary conditions on the components of R. 



An evolution model for (3.4) can be obtained by substituting (3.5) into (3.6) which 
yields 



(f){K,,y,t) = A{k) (f){K,y,t) + B{k,) d{n,y,t) , 
v{K,,y,t) = C{k,) (f){K,,y,t) + D(a^) d{K,,y,t), 
where the operators A, B, C, and D are defined in Appendix [Aj 



(3.7a) 
(3.76) 



3.3. Spatio-temporal frequency responses 

Application of the temporal Fourier transform yields the frequency response operator for 
system (3.7) 

H(A^,cj;/3,We,L) = C (a^) (icj I - A (a^))"^ B (a^) + D (a^) , (3.8) 

where u is the temporal frequency, and I is the identity operator. For a stable sys- 
tem (3.7), (3.8) can be used to characterize the steady-state response to harmonic input 
signals across spatial wavenumbers k, and temporal frequency uj. Namely, if the input d 
is harmonic in x, and t, i.e., 

d{x,y,z,t) = d(y)e^(^»-+^^^ + ^*), 

with d{y) denoting some spatial distribution in the wall-normal direction, then the output 
V is also harmonic in x, 2:, and t with the same frequencies but with a modified amplitude 
and phase 

vix,y,z,t) = ([H(fc,,fc„w)d](y))e'(^^^ + ^^^ + ^*), 
where the amplitude and phase are precisely determined by the frequency response at the 
input frequencies (kx^ kz^oj). Note that we have dropped the dependence of the frequency 
response operator on We, /3, and L for notational convenience. 

The nth singular value of the frequency response operator H is determined by 

alin) = A„(H*H), 
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where An(-) denotes the nth eigenvalue of a given self-adjoint operator and H"^ is the 
adjoint of H. For any {kx^ k^^oj)^ crmax(H) = max^ crn(H) determines the largest amplifi- 
cation from d to V. Furthermore, the temporal supremum of the maximal singular value 
of H determines the Hoo norm of system ( |Zhou et a/.||1996D 



G{K;P,We,L) 



sup cr^ax (H {f^^ ^; We, L)) . 



This measure of input-output amplification has several appealing interpretations: 

(a) for any {kx^kz), the Hoo norm represents the worst-case amplification of purely 
harmonic (in and t) deterministic (in y) disturbances. This worst-case input-output 
gain is obtained by maximizing over input temporal frequencies (sup over u) and wall- 
normal shapes (maximal singular value of H); 

(b) in the temporal domain, the Hoo norm represents the energy gain from forcing to 
velocity fluctuations 



sup / \ : 

Edin) < 1 ^d('^) 



where E^{k) denotes the A^-parameterized energy of velocity fluctuations, i.e., 

pOO pi 

E^{n) = J J w*{n,y,t)w{K,y,t)dydt. 

In other words, for a unit-energy forcing, G{k) captures the largest possible energy of 
velocity fluctuations across wavenumbers and 

(c) at any {kx^ kz), the inverse of the Hoo norm quantifies the size of an additive un- 



structured modeling uncertainty F that can destabilize generator A in (3.7). As described 
in § [2j large H^o norm indicates small stability margins (i.e., low robustness to model- 
ing imperfections). For systems with poor robustness properties, even small modeling 
uncertainties (captured by operator F) can lead to instability of operator A + F. 

We also note that the frequency response of system (3.7) can be further decomposed 
into 3x3 block-operator form 



u 




V 




w 





Hnl H^2 H^3 

H^i H^2 H^3 

H^(;l H^2 H^3 





di 








_ ds _ 



(3.9) 



This form is suitable for identifying forcing components that introduce the largest am- 
plification of velocity fluctuations. In (3.9), H^j maps dj to s, and 



Gsj{K,]f3,We,L) = sup cr^^^ (H^^- (k,cj;/3, We,L)) , s = {u,v,w}, j = {1,2,3}. 

UJ 

The finite-dimensional approximations of the underlying operators are obtained us- 
ing Matlab Differentiation Matrix Suite CWeideman & Reddy 2000), which utilizes 



pseudospectral methods to approximate differential operators. After discretization in the 



wall- normal direction, each component in (3.9) becomes an x matrix, where N de 



notes the number of Chebyshev collocation points in y. All computations are performed 
in Matlab and grid-point convergence is confirmed by running additional computations 
with larger number of grid points in y. 

After discretization in y, the H^o norm of the frequency response matrix can, in prin- 
ciple, be computed by determining crmax(H(cj)) for many values of u and by choosing 
the resulting maximum value. However, there are two obvious problems associated with 
such a method: difficulty in determining the range and spacing of the temporal frequen- 
cies, and the large number of computations. To avoid these issues, [Boyd, Balakrishnan &[ 
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Kabamba (1989) devised a bisection method that can efficiently compute the Hoo norm. 



Furthermore, Bruinsma & Steinbuch ( 1990 ) introduced a fast algorithm that utilizes an 
efficient method of choosing the temporal frequency for computing the Hoo norm. This 
fast algorithm is utilized in our computations and it is based on the relation between 
the singular values of the frequency response matrix and the eigenvalues of a related 
Hamiltonian matrix. 

All of our results are conffimed by additional frequency response computations that 
utilize the integral formulation of (3.7). This is accomplished by rewriting the evolution 
equations (3.7) into an equivalent two-point boundary value problem and then refor- 
mulating it into a system of integral equations. The procedure for achieving this along 
with easy-to-use Matlab source codes is provided in Lieu fc Jovanovic^(2Q ll). This new 
paradigm for computing frequency responses utilizes the Chebfun computing environ- 
ment ( Trefethen et al [2011 ) and it exhibits superior numerical accuracy compared to 
conventional numerical schemes. 



4. Frequency responses of 3D velocity fluctuations 

In this section, we study the dynamics of three-dimensional velocity fluctuations in an 
inertialess shear-driven channel flow. In particular, we examine the worst-case amplifl- 
cation of deterministic disturbances and identify the corresponding wavenumbers that 
contain most energy. Our analysis shows that velocity fluctuations with large stream- 
wise and 0{1) spanwise length scales display the highest sensitivity to disturbances (and 
consequently the lowest robustness to modeling imperfections). We further utilize the 
component- wise frequency responses ( Jovanovic fc Bamieh|2005| to identify forcing com- 
ponents that have the strongest influence on the velocity fluctuations. In strongly elastic 
flows, we demonstrate that the wall- normal and spanwise forces have the highest impact, 
and that the streamwise velocity is most amplifled by the system's dynamics. 

In the remainder of this section, we set the viscosity ratio to /3 = 0.5 and study the 
effect of the Weissenberg number, We, and the maximum dumbbell extensibility, L, on the 
kinetic energy density. Figure [4] shows the worst-case ampliflcation of velocity fluctuations 
triggered by ah three body forces in flows with We = {10, 50, 100} and L = {10, 50, 100}. 
For L = 10, the spatial frequency responses display low-pass fllter features commonly 
seen in flows governed by viscous dissipation, with the peak ampliflcation taking place 
at low wavenumbers. Furthermore, this spatial distribution remains almost unchanged 
as We increases from 10 to 100 (cf. flgures [4(a)| 4(d)[ and |4(g)[ ). However, as L increases 
to 50, the velocity fluctuations become more amplifled with Weissenberg number and a 
dominant peak starts to appear in an isolated region around kx ^ O(10~^) and ~ 
0(1); see flgures 4(e)| and |4(h)[ As L increases to 100, ampliflcation with We increases 
even further. This indicates that in flows with large enough L and We, streamwise- 
elongated and spanwise-periodic flow fluctuations are the most amplifled by deterministic 
body forces. Thus, in strongly elastic channel flows without inertia, streamwise-constant 
and nearly streamwise-constant fluctuations with a preferential spanwise length scale are 
most sensitive to external disturbances. 

We next study the component-wise frequency responses that quantify ampliflcation 
from different forcing to different velocity components. This analysis facilitates identifl- 
cation of forcing components that are most effective in amplifying velocity fluctuations. 
In Couette flow with We = 50 and L = 50, flgure [5] shows the worst-case ampliflcation of 
the 9 frequency response components in (3.9). We see that the streamwise and spanwise 
velocity fluctuations are the most amplifled. Furthermore, the maximum ampliflcation of 
the streamwise velocity fluctuations triggered by d2 and occurs around k^, ^ O{10-^) 
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Figure 4. Worst-case amplification from d to v in Couette flow with = 0.5, 
We = {10, 50, 100}, and L = {10, 50, 100}: first row. We ^ 10; second row. We = 50; and third 

row. We = 100. 



and kz « 0(1), respectively. This illustrates that streamwise velocit y fluc tuations are 
responsible for the most amplified region {kx ~ 10~^, kz ~ 2) in figure 4(e) In contrast, 
the square region around kx ^ and kz ^ in figure 4(e)| arises from the responses of 
streamwise and spanwise velocity fluctuations to di and ds, respectively. We note that the 
wall-normal velocity experiences negligible amplification compared to that of the other 
two velocity components and hence, does not contribute to the large energy amplification 
in inert ialess flows. 

The above results clearly illustrate the dominance of streamwise-constant and nearly 
streamwise-constant velocity fluctuations in strongly elastic Couette flow without iner- 
tia. The streamwise velocity is most amplified by disturbances, and this large response 
is caused by the wall-normal and spanwise body forces. Even though this section pro- 
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Figure 5. Component- wise worst-case amplification from dj to s in Couette flow with 
s = {u,v,w}, j = {1,2,3}, /3 = 0.5, We = 50, and L = 50. The symbol (•) identifies the 
largest value of the corresponding plot. 



vides useful insight into the dynamics of inertialess channel flows of viscoelastic fluids, 
the scaling of energy amplification with We and L cannot be deducted from our compu- 
tations. For streamwise-constant velocity fluctuations, this issue is addressed in § |5] by 
developing explicit analytical expressions that quantify the dependence of the worst-case 
amplification on We, L, /3, and kz- 

5. Dynamics of streamwise-constant velocity fluctuations 

Motivated by the observation that in strongly elastic flows the most amplified velocity 
fluctuations have large streamwise length-scales, we next examine the linearized model 
for fluctuations without streamwise variations, i.e., at kx = 0. We use this model to 
establish an explicit scaling of the components of the frequency response operator with 
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We and L, and to explain the observations made in § [4j Since the largest amplification 
takes place at low temporal frequencies, we use analytical developments to show that in 
flows without temporal variations the worst-case amplification from d2 and ds to u scales 
linearly with the first normal stress difference Ni of the nominal flow. Consequently, 
this worst-case amplification scales quadratically with We as L ^ oo and quadratically 
with L as We oo. Therefore, even in flows with infinitely large polymer relaxation 
times finite extensibility of polymer molecules limits the largest achievable amplification. 
Furthermore, the worst-case amplification from all other forcing to all other velocity 
components is both We- and L-independent. We also present the spatial structures of 
the forcing and velocity fluctuation components that contribute to the above mentioned 
unfavorable scaling with Ni , and demonstrate that the key physical mechanism involves 
interactions of polymer stress fluctuations in the z)-plane with base shear. 



5.1. Dependence of worst-case amplification on We and L 

For fluctuations without streamwise variations, i.e. at kx = 0, equations (3.6) and (3.5) 
simplify to 



01 

02 



-Fii 01 

-F22 02 

1 - (3 



F2I 01 



F2v V + F2r, V, 



1 - (3 



C^2 02 



u 




V 




w 





C„jj 

1 



1 



D^i di 



1 



(5.1a) 
(5.16) 

(5.1c) 
{b.ld) 

(5.1e) 



where (f)i = [ r22 ^23 ^33 ]^ ^ 4^2 — [ ^13 ^12 ^11 ]^ 7 is the wall-normal velocity, 
and T] is the wall-normal vorticity. On the other hand, the operators in (5.1 ) are given by 
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An evolution representation of the streamwise-constant model (5.1) can be obtained 



by eliminating the components of the conformation tensor from the equations. This is 



achieved by substituting the temporal Fourier transforms of (5.1a) and (5.16) into (5.1c 



and (b.ld) and taking the inverse temporal Fourier transform of the resulting equations. 
We will show that this representation leads to convenient analytical expressions for the 
frequency response operator. 



For streamwise-constant fluctuations the frequency response operator H in (3.9) sim- 
plifies to 



(5.2) 



u 




H^i 


Hn2 


Hn3 




di 
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H^3 
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w 







tiw2 


H^(;3 




ds 



where the operators Hgj are given by 
ll,2{k,,co;l3,We,L) 

H„3(fc.,w;/3,We,L) 
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iuj -i 


-f 


iujP - 


f / 
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12 (iuj(3 + f) (Co - + ic^Ci) 
Wef{/3-l) 2We^{/3-l) 3iujf-uj'^ 



^2 



/ 



£2 
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= 2/ 



{icv/3 + f) (Co - w2 + iwCi) 
2We 



(5.3) 
(5.4a) 

A-'dyy, (5.46) 



L2 



(5.4c) 



In (5.3), we have successfully separated the temporal and spatial parts of the frequency 
response operators from d2 and ^3 to v and w. The absence of inertia induces simple tem- 
poral dependence of Hgj with {s = v, w; j = 2, 3} and facilitates analytical determination 
of the temporal frequency uj at which the largest worst-case amplification takes place. 
These four frequency response operators exhibit high-pass temporal characteristics, with 
the peak amplification taking place at infinite frequency, w = 00. The addition of a small 
amount of inertia would introduce roll-off at high temporal frequencies, thereby shifting 
the peak amplification to finite temporal frequency ( Jovanovic fc Kumar||2011 ). 

These observations allow us to obtain explicit expressions for the worst-case amplifi- 
cation from d2 and ds to v and w. For example, the worst-case amplification from d2 to 
V is given by 

G„2(fc.;/3) 



sup cr^ax {tlv2{kz,uj;We,l3,L)) 



(I//32) <tLx(D„2) 

(1//3') 5.2(fc.), 
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Figure 6. Functions characterizing worst-case amplification from d2 and to v and w at 
kx — 0; gsj{kz) with {s — v,w] j — 2,3}. 



where the function ^^2, which is independent of We, /3, and L, captures the spanwise 
frequency response (from d2 to v). A similar procedure yields the following expressions 
for the worst-case amplification from d2 and ds to v and w 



(5.5) 



where the functions g represent the We-, and L-independent spanwise frequency 
responses. We note that, in Couette flow without inertia, the worst-case amplification of 



the four components in (5.5) is equivalent for Oldroyd-B and FENE-CR fluids. 

The functions gsj{kz) with {s = w] j = 2, 3}, are shown in figure |6] We see that 
and gys exhibit similar trends with peaks at kz ~ we also note that gw2 = gvs- In 

contrast, ^^3 has a low-pass shape with maximum occurring at kz =0. The peak value 
of this function is about four times larger than the peak values of ^^2 and ^^3. 

We next examine responses of streamwise velocity to different forcing components. We 
first analyze the temporal characteristics of the frequency response operators Huj with 
j = {1, 2, 3}. Since the largest amplification of u arising from d takes place at kz ~ 
in figure [t] we show the temporal frequency dependence of cr^ax(iiuj) foi" kz = 1.5, 



/3 = 0.5 and L = 10. Figure |7(a)| shows high-pass temporal features of (Jmax (H^^i) with 
its maximum value taking place at = 00; furthermore, this peak value is independent 
of the Weissenberg number and the maximum extensibility of the polymer molecules L. 
In contrast, (Jmax (H^2) and (Jmax (H^3) have low-pass characteristics and attain their 
largest values (which depend on both We and L) at low temporal frequencies. 

Figure [8] shows the /C;2-dependence of the functions Guj that quantify the worst-case am- 
plification from different forcing components to the streamwise velocity for L = {10, 100} 
and for multiple values of the Weissenberg number. Compared to the wall- normal and 
spanwise velocity fluctuations, streamwise velocity is more amplified by disturbances; 
cf. figures |6] and [s] As evident from figures |8(a)| and |8(d)[ Gui has high values for low 
spanwise wavenumbers and is independent of both We and L. On the other hand, Gu2 
and Gu3 achieve their peaks at kz ~ 0{1). For a fixed value of L these two functions 
increase with the Weissenberg number, and the largest amplification takes place in the 
limit of infinitely large We. In contrast to the Oldroyd-B fluids, finite extensibility of the 
nonlinear springs in the FENE-CR model induces finite values of Gu2 and Gus even at 
arbitrarily large Weissenberg numbers. Furthermore, for a fixed value of We, worst-case 
amplification increases with L and the largest amplification is obtained in the Oldroyd-B 
limit (as L 00). Analytical explanation for these observations is provided below. 




10° 10^ lO'' 10"^ 10° 10^ 10"^ 10° 10^ 
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Figure 7. Maximum singular values of the frequency response operators from dj to as a 
function of uj in streamwise-constant Couette flow with j = {1,2, 3}, kz = 1.5, /3 = 0.5, L = 10, 
and We — [10, 100]. The symbol (o) identifies the peak values of the corresponding curves. 
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Figure 8. Worst-case amplification from dj to u in streamwise-constant Couette flow: first row, 
L = 10 and We = [10,100]; and second row, L = 100 and We = [10,1000]. The symbol (•) 
shows the worst-case amplification from dj to u in the limit of infinitely large We. 



We next present explicit expressions for the worst-case amplification from different 
forcing components to the streamwise velocity. We summarize our major findings here 
and relegate derivations to Appendix |B] We first consider the worst-case response of the 
streamwise velocity in the presence of the streamwise body forcing. For all values of We 
and L, our computations indicate that the worst-case amplification from di to u takes 
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Figure 9. Functions characterizing worst-case amplification from d2 and ds to u at kx 

oj = 0: gu2 (o); and gu3 (solid). 



and 



place at = oo. Consequently, in the limit of infinitely large uj we have 



G„i(fc,;/3) 



lim 



.{■H^i{k,,uj;p,We,L)) 



(5.6) 



where gui is the spanwise frequency response from di to u. Note that Gui is independent 
of both We and L which is in agreement with the observations made in figures 8(a)| 
and [8(d)| 

Derivation of the analytical expressions for Gu2 and Gus is more challenging because 
the worst-case amplification of u arising from ^2 and ds depends on both We and L. 
However, since our computations demonstrate that the worst-case amplification from 
d2 and ds to u takes place at low temporal frequencies, the essential features can be 
captured by analyzing the corresponding frequency responses at cj = 0. Thus, the worst- 
case amplification of u caused by d2 and ds can be reliably approximated by 



Guj{k,;f3,We,L) ^ a, 



(H,,(fc„0;/3,iye,L)) 
{We/f)'{l 
1 



= 0^1(1 



P) guj{k^) 

uj{kz)^ j = {2,3}, 



(5.7) 



where the functions guj with j = {2,3} quantify the spanwise frequency responses from 
d2 and ds to u; see figure |9] Equation ( |5.7[ ) shows that the worst-case amplification of the 
streamwise velocity fluctuations scales linearly with the first normal stress difference Ni 
of the nominal Couette flow. Hence, even in the absence of inertia, velocity fluctuations 
can experience large amplification in flows with large normal stress difference. 

Furthermore, table[l]shows explicit expressions for Gu2 and GuS in the limit of infinitely 
large We (or infinitely large L). In these two cases, the first normal stress difference of 
the base Couette flow is given by 



lim Ni = 2We^ 



lim A^i 



We- 



For Oldroyd-B fluids polymer molecules are modeled by infinitely extensible linear springs 
and both Gu2 and Gus scale quadratically with the Weissenberg number. On the other 
hand, as We — ^ 00, Gu2 and Gus scale quadratically with the maximum extensibility 
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L — > oo I We — > oo 



Gu2 (/c.; /3, •, •) I We^ (1 - gu2{k,) I 0.5 (1 - gu2{k,) 



Gu3 (/c.; /3, •, •) I We^ (1 - /Sf gu3{k^) | 0.5 (1 - gu3{k^) 

Table 1 . Worst-case amplification of streamwise velocity fluctuations arising from the wall-nor- 
mal and spanwise forces in the limit of infinitely large maximum extensibility of polymer chains 
L or infinitely large Weissenberg number We. 



of the polymer molecules L. We conclude that - even for infinitely large polymer relax- 
ation times - energy amplification of velocity fluctuations in inertialess Couette flow of 
viscoelastic fluids is bounded by the maximum extensibility of nonlinear dumbbells. 



5.2. Dominant flow structures 

In this section, we present the spatial structures of the wall-normal and spanwise body 
forces that induce the largest amplification in streamwise velocity. We also discuss the 
z)-dependence of the resulting streamwise velocity fluctuations. These structures are 
purely harmonic in the spanwise direction with period determined by the value of kz at 
which the functions gu2 and gus attain their maxima {kz ~ 2.5 and kz ~ 1.5, respectively). 
Furthermore, since the worst-case amplification occurs at = 0, these structures are 
constant in time and their wall-normal profiles are determined by the principal singular 
functions of the frequency response operators that map d2 and ds to u. 

Fluctuations in d2 and ds that lead to the largest amplification of u are shown in 
figures p"0( a) and 10(b) The wall- normal forcing is symmetric with respect to the channel 



centerline with the peak value located at the center of the channel; in contrast, the 
spanwise forcing is antisymmetric with respect to the channel centerline. 

Figures p!Q(c) and 1 10 (d)] illustrate the spatial structures of streamwise velocity induced 
by the body forcing fluctuations shown in figures pi) ( a) | and [l (b ) | Both body forces yield 
a symmetric response in u with vortices occupying the entire channel width. We note that 
the dominant flow structures shown in figure 10 do not exhibit significant deviation with 
/3, We, and L. Furthermore, we observe striking similarity between these flow structures 
and flow structures resulting from the analysis of stochastically forced Couette flow of 
Oldroyd-B fluids ( Jovanovic fc Kumar||201 1 ) . Although these spatial structures may not 
match the full complexity of flow patterns produced in experiments and direct numerical 
simulations, our analysis identifies dynamical features that are likely to play an important 
role in shear-driven channel flows of viscoelastic fluids. 



5.3. Physical mechanisms 

We next discuss the physical mechanisms responsible for strong influence of wall-normal 
and spanwise forces on streamwise velocity fluctuations in inertialess channel flows. As 
shown in Appendix [Bj in the absence of streamwise forcing and streamwise variations in 
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Figure 10. First row: body forcing fluctuations in d2 and ds leading to the largest amplification 
of u. Second row: streamwise velocity fluctuations obtained by forcing the inert ialess Couette 
flow with the body force fluctuations shown in (a) and (b), respectively. 



flow fluctuations, u{y^ t) evolves according to 
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Therefore, even in the absence of inertia, the source term in the evolution equation for 
streamwise velocity is provided by the interactions between polymer stress fluctuations 
(in the wall-normal/spanwise plane) with the base shear U' . Furthermore, the finite ex- 
tensibility of the polymer chains introduces additional source terms that are not present 
in the Oldroyd-B model; physically, these additional terms originate from: (i) the in- 
teraction between the streamwise shear component of the conformation tensor ri2 with 
the base shear; and (ii) the wall-normal gradient of the streamwise component of the 
conformation tensor rn. Relative to the terms that are already present in the Oldroyd-B 
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model, the influence of these additional terms (that arise from temporal changes in rn) 
is much weaker. 



We have demonstrated in § 5.1 that the the essential features of the worst-case ampli- 
fication from d2 and ds to u can be captured by analyzing the corresponding frequency 



responses at a; = 0. At zero temporal frequency, equation (5.8) simplifies to a static- 



in-time relation between the streamwise velocity and the fiuctuating components of the 
polymer stress tensor. 



Au 



0^ (1 - (3) {dy {U't22) + (U't2^) ; 



(5.9) 



Furthermore, T22 and r23 are proportional to the spatial gradients of the (?/, z)-plane 
streamfunction (i.e., v = dz^l^^ w = —dytjj)^ 



T22 = '^dyz'il), r23 = - {dyy - dzz)"^, 
and iIj is induced by the action of the wall-normal and spanwise forces 

d2 



^ = [ dy ] 



Finally, by substituting (5.10) into (5.9) we obtain the following expression 



Au = 



0Vi/2 (i-/?)a,AV' 

7^(1-/3)5,^,, 



(5.10) 



(5.11) 



(5.12) 



that relates fiuctuations in the streamwise velocity u and the streamwise vorticity ujx in 
inertialess Couette fiow of FENE-CR fiuids without streamwise and temporal variations 
(i.e., at /Ca^ = and a; = 0). 

This demonstrates that 0{1) fluctuations in streamwise vorticity induce 0{\fNi) fluc- 
tuations in streamwise velocity through a viscoelastic equivalent of the well-known lift-up 
mechanism. In contrast to Newtonian fluids, where vortex tilting induces large amplifl- 
cation, the lift-up mechanism in viscoelastic fluids originates from interactions between 
polymer stress fluctuations in the (?/, z)-plane with background shear (jJovanovic fc Ku- 



mar 



2011). In the absence of inertia, a static-in-time momentum equation relates the 



wall-normal and spanwise velocity fluctuations (and consequently the streamwise vortic- 
ity) to the polymer stress fluctuations r22, T23, and raa. Interactions of these polymer 
stress fluctuations with background shear induce the energy transfer from the mean flow 
to fluctuations and redistribute momentum in the (?/, z)-plane through a movement of 
the low speed fluid (away from the wall) and the high speed fluid (towards the wall). 
This momentum exchange is responsible for the generation of alternating regions of high 
and low streamwise velocity (relative to the mean flow), and it is facilitated by large 
normal stress difference TVi, low viscosity ratios /3, strong base shear U' ^ and strong 
spatial variations in streamwise vorticity fluctuations. As in streamwise-constant inertial 
flows of Newtonian fluids, this ampliflcation disappears either in the absence of spanwise 
variations in flow fluctuations or in the absence of the background shear. 



6. Dynamics of streamwise-constant polymer stress fluctuations 

Although we have so far conflned our attention to the dynamics of velocity fluctu- 
ations, it is worth noting that polymer stress fluctuations can also experience signifl- 
cant ampliflcation even in the absence of inertia. Since our computations (not shown 
here) demonstrate that largest responses in polymer stress fluctuations are induced by 
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streamwise- const ant deterministic forcing, we next examine the responses from body 
forcing to polymer stress fluctuations without streamwise- variations. We use analytical 
developments to show that the wall-normal and spanwise forces induce the largest am- 
plification of the polymer stress fiuctuations. The worst-case amplification obtained in 
the presence of these two forcing components takes place at a; = and it is proportional 
to: Ni for tis and ri2; and TVf / (l + 7Vi/Z^) for rn. Furthermore, the worst-case am- 
plification from all forcing components to r22, T23, and rsa is independent of /3, We, and 
L. We also illustrate that the worst-case amplification from d2 and ds to rn scales 

• quartically with We as L ^ 00; 

• quartically with L as We 00. 

Following a sequence of straightforward algebraic manipulations, the frequency re- 
sponse operator G that maps body forcing fiuctuations di, ^2, and ds to polymer stress 
fiuctuations can be expressed as 



T22 




" 


G12 


Gi3 


T23 







G22 


G23 


T33 







G32 


G33 


ri3 




G41 


G42 


G43 


ri2 




G51 


G52 


G53 


^11 






^62 


^63 



di 



(6.1) 



where the streamwise- const ant frequency response operators G^j are given in Appendix[C| 
We note that all components of the frequency response operator in (6.1) exhibit roll- 
off at high temporal frequencies, thereby indicating that the largest singular value of 
each component of G peaks at finite temporal frequency. In particular, the worst-case 
amplification from d2 and ds to r22, T23, and r^s takes place at cj = and is determined 
by the following We-, and L-independent functions 



Gi2{h) 

G22{h) 

Gs2{h) 



G23{h) 

Gssih) 



922{k^) g23{kz) 



(6.2) 



Figure pT] shows the functions gij with {£ = 1, 2, 3; j = 2, 3} that quantify the spanwise 
wavenumber dependence of the respective frequency response operators. We note that 
932 = 9i2 and gss = gis- From figure pT(a)[ we see that gi2 and gis decay to zero at 
both low and high wavenumbers with the maximum values occurring at kz ~ 2.3 and 
kz ~ 1.5, respectively. Similarly, function ^22 decays to zero at both low and high values 
of kz and it achieves two peaks at kz ^ lA and kz ^ 4.0. On the other hand, ^23 displays 
low-pass behavior and the maximum value of this frequency response is approximately 
three times larger than the maximum values of other responses in figure [TT] 

We next analyze the temporal frequency responses of the operators G^j with {£ = 
4, 5,6; j = 1,2,3} that map different forcing components to fiuctuations in ri3, ri2, and 
Til. Figure 12 shows a;-dependence of (Jmax (G^j) for kz = 1.5, (3 = 0.5, and L = 10. 
Figures 12(a) - 12(c) show that cTmax (^4^) with j = {1,2,3} achieve their respective 
peaks at = 0. Furthermore, while the peak value of (Jmax (G41) does not depend on the 
Weissenberg number, the peak values of (Jmax (G42) and (Jmax (G43) increase with We. 

The frequency responses G^j with {£ = 5,6;j = 2,3} that quantify amplification 
from d2 and ds to ri2 and rn exhibit similar low-pass characteristics. On the other 
hand, figures 12(d) and 12(g) show that (Jmax (G51) and (Jmax (Gei) achieve their peak 
values at non-zero temporal frequencies and that these values increase as We increases. 
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9i2{kz), 9i3{kz)^ 932{kz), gssikz) ^22 (^z), g23{kz) 




We see that the forcing components in the wah-normal and spanwise directions induce 
larger amphfication of polymer stress fluctuations compared to the streamwise forcing. 
Furthermore, the streamwise component of the polymer stress tensor rn experiences the 
largest amplification. 

Following a series of algebraic manipulations, it can be shown that the worst-case 
amplification from ^2 and ds to ri3, ri2, and rn takes place at = and is given by 



G^j{h;P,We,L) = {2 + pf g^,{k,), 

Gej{k,;l3,We,L) = (1 + gsife), j = {2,3}. 

(1 + Ni/L^) 



(6.3) 



The functions gij with {^ = 4, 5, 6; j = 2, 3} represent the We- and L-independent span- 
wise frequency responses from the wall-normal and spanwise forces to ri3, ri2, and rn. 



Equation (6.3) shows that the worst-case amplification of ris and ri2 is proportional to 

v2 



Ni. On the other hand, the worst-case amplification of rn scales as Nf/ (l -\- Ni/L' 



'2\ 



Figure 13 shows the A:^ -dependence of the functions g£j for {£ = 4,5,6; j = 2,3}. We 
note that ^52 = g62 and ^53 = g^s. Function ^43 has low-pass shape and it peaks at 
kz = 0. We also notice band-pass features of ^52, ^53, and ^42, with the peak values 
occurring at kz ~ 2.4, kz ~ 1.6, and kz ~ 1.4, respectively. 

Table [2] summarizes the worst-case amplification of ri3, ri2, and rn arising from the 
wall- normal and spanwise forces in the limit of infinitely large L (or infinitely large We) . 
For Oldroyd-B fluids (i.e., as L ^ 00), both G^j and G^j scale quadratically with the 
Weissenberg number, and Gqj scales quarticly with We. On the other hand, as We 00, 
both G^j and G^j scale quadratically with L and Gqj scales quarticly with L. This 
demonstrates profound influence of ^2 and ds on rn in strongly elastic shear flows of 
viscoelastic fluids. Thus, even in the absence of inertia viscoelastic shear flows with 
large polymer relaxation times and large extensibility of polymer molecules exhibit high 
sensitivity to disturbances and low robustness to modeling imperfections. 
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(jj UJ 



(g) (h) (i) 

Figure 12. Maximum singular values of the streamwise-constant frequency responses operator 
from dj to ns (G4j), T12 (Gsj), and m (Gej) as a function of uj for j = {2,3}, kz = 1.5, 
/3 = 0.5, L = 10 and We = [10, 100]. 

7. Concluding remarks 

In this study, we have examined non-modal amplification of disturbances in inertialess 
Couette flow of viscoelastic fluids using the FENE-CR model. The amplification is quanti- 
fied by the maximal singular values of the frequency response operators that map sources 
of excitations (body forces) to the quantities of interest (velocity and polymer stress fluc- 
tuations) . Spatio-temporal body forcing fluctuations are assumed to be purely harmonic 
in the horizontal directions and time, and deterministic in the wall-normal direction. Our 
three-dimensional component-wise frequency response analysis of the FENE-CR model 
sets the current paper apart from prior works which study the transient growth of veloc- 
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^42 (/^^z; 0.5), gA3{kz;0'5) 952{kz), g^si^z)^ 9Q2{kz)^ gesi^z) 




j = {2,3} I L — ^oo I We — ^ oo 

G4j{k.;Prr) I We^g4j{k.;P) I 0.5 g4j{k.; P) 



Gsj {k.;l3, •, •) I We^ (2 + /Sf g,j (k.) I 0.5 (2 + fif g,j (k.) 



G6j{k,;l3,-,-) I 4We^{l + 2l3f g6j{k,) | 0.25 (1 + 2/3)' ^e, (/c.) 

Table 2. Worst-case amplification of ns, ri2, and m arising from (i2 and in tfie limit of 
infinitely large maximum extensibility or infinitely large Weissenberg number. 



ity and polymer stress fluctuations in inertialess flows (Jovanovic & Kumar 2010) and 
non-modal amplification of stochastic disturbances in elasticity-dominated flows with 
non-zero inertia ( Hoda et a/.||2008 , 2009; Jovanovic fc Kuniar||2011 ) using the Oldroyd- 
B model. We have shown that streamwise-elongated flow structures are most amplified 
by disturbances. Furthermore, the component-wise frequency responses reveal that the 
wall-normal and spanwise forces have the strongest impact on the flow fluctuations, and 
that the influence of these forces is largest on streamwise components of velocity and 
polymer stress fluctuations. 

For streamwise-constant fluctuations, we have established analytically that the largest 
amplification of the streamwise velocity and streamwise component of the polymer stress 
tensor is proportional to the first normal stress difference of the nominal flow Ni and to 
Ni/ (l + 7Vi/Z^) , respectively. This largest amplification is caused by wall- normal and 
spanwise forcing fluctuations and it takes place at low temporal frequencies and 0{1) 
spanwise wavenumbers. Using our analytical developments we have also shown that this 
worst-case amplification of u and rn respectively scales as (i) O(We^) and (9(We^) in 
the Oldroyd-B limit (i.e., as L ^ oo); and (ii) 0{L'^) and 0{L^) in the limit of infinitely 
large Weissenberg number. We thus conclude that in the presence of large polymer re- 
laxation times and large extensibility of the polymer molecules, the velocity and polymer 
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stress fluctuations can experience significant amplification even when inertial effects are 
completely absent. The underlying physical mechanism involves interactions of polymer 
stress fluctuations with a background shear, which induces a viscoelastic analog of the 
vortex tilting mechanism that is responsible for large amplification in inertial flows of 
Newtonian fluids. 

It is worth noting that, in the limit of infinitely large We, the worst-case amplification 
of both velocity and polymer stress fluctuations is bounded by the maximum extensibility 
of the polymer molecules. This is in contrast to Oldroyd-B fluids where infinite extensi- 



bility allows the amplification of disturbances to grow unboundedly with We (Jovanovic 
fcKumar|2QlT] ). Our new observations demonstrate that high sensitivity to disturbances 



and low robustness to modeling imperfections are reduced by finite extensibility of non- 
linear dumbbells. Thus, both large polymer relaxation times and large extensibility of 
the polymer molecules are needed to achieve large amplification of velocity and polymer 
stress fluctuations in inertialess channel flows of viscoelastic fluids. 

The present work extends recent efforts ( Hoda et a/.||2QQ8 , 2009 Jovanovic & Kumar 



2010, 2011) that examine possible mechanisms for triggering transition to elastic turbu- 



lence in channel flows of viscoelastic fluids. In addition to providing insight into worst- 
case amplification of velocity and polymer stress fluctuations in inertialess flows, we also 
demonstrate the importance of uncertainty quantification in flows of viscoelastic fluids. 
Our analysis shows high sensitivity of inertialess flows of viscoelastic fluids to external 
disturbances. Unfavorable scaling of the worst-case amplification of flow fluctuations with 
We and L indicates that small-in-norm modeling imperfections can destabilize nominally 
stable flows. Hence, stability margins of inertialess channel flows of viscoelastic fluids 
decrease significantly with an increase in the Weissenberg number and the maximum 
extensibility of the polymer chains. This uncertainty may arise from inevitable imperfec- 
tions in the laboratory environment or from the approximate nature of the constitutive 
equations. Our observations regarding model robustness also have important implica- 
tions for numerical simulations, where numerical and/or roundoff errors may cause the 
simulated dynamics to differ from the actual dynamics. 

The present findings suggest a plausible mechanism for transition to elastic turbulence 
in channel flows of viscoelastic fluids. Large amplification of disturbances induces forma- 
tion of streamwise streaks whose growth can put the flow into a regime where nonlinear 
interactions are no longer negligible. These nonlinear interactions can then induce sec- 
ondary amplification ( Schoppa fc Hussain||2002 ) or secondary instability ( Waleffe||1997 ) 
of streamwise streaks, their breakdown, and transition to a time-dependent disordered 
flow and elastic turbulence. To understand possible routes for the transition to elastic 
turbulence, it is essential to track later stages of disturbance development by considering 
nonlinearities in the constitutive equations and their interplay with streak development 
and high flow sensitivity. Our ongoing efforts are directed toward examining sensitivity 
of the streaks to three-dimensional disturbances. We also intend to study the presence 
of a self-sustaining mechanism (proposed for Newtonian fluids by Waleffe] ( 1997[ )) and 
to numerically track later stages of disturbance development in strongly elastic channel 
flows of viscoelastic fluids. 
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Appendix A. The underlying operators in 3D shear-driven channel 
flow of FENE-CR fluids 



In this appendix, we define the underlying operators appearing in (3.5)-(3.7) for a 
shear-driven channel flow. Operators {C,D} in (3.5) are given by 

(1-/ 



— 



13 



(1 



1 

1 



C^i Cy2 ] , = ^ [ D^i D^2 ] , 



D^l D^2 DryS J , 



where 

C^i = 

C^i = 

^vl,l = 

^vl,3 = 

^v2,2 = 

Cryl,! = 

Cr7l,3 = 

^v2,2 = 

D,i = 

D^i = 



A ^ [ C^i^i C^i^2 ^vl,3 ] , C^2 — A ^ [ C^2,l Cv2,2 Cv2,3 ] , 
A"-*- [ C^i^i C^i^2 C^i^3 ] , C^2 = A"-*- [ C^2,l C^2,2 C^2,3 ] 5 



f klf 2Wekl 
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fc^), c 



't;2,l 



2 kzkx f 
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~W^ 



(dyy + /C^ 



9„ H —[Oy 
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(2T^e/L2)fc,A:, - {ikJ/D)dy, C^i,2 = (ifc, //W^e) 
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We, 



t^j/) C„2,3 — 



2 We 



We 



ik^A-^dy, D„2 = PA-^, D„3 = i/c^A-^Sj,, 
-ik,A-\ D^2 = 0, D^3 = ik^A-\ 



Here, fc^ = ^2 _|_ ^2^ j _ /\ = ^j^^ _ ^j^j^ Dirichlet boundary conditions, 



yyyy 



2k'^dyy + k^ with both Dirichlet and Neumann boundary conditions. 



The F-operators appearing in (3.6) are determined by 
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2e {dy + ik^Ri2) 



ikz {dy 



We kx kz Ri2 dy 
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\kxdy - klRi2 



7 -^277 



We 
1^ 



kl + \kx Ri2 dy - kl Rii 

kx kz Ri2 - '^kz dy 
2 {kx kz Rii - ikz R12 dy) 



2 kx kz 

ikz R12 dyy - kxkz (1 + ^11) dy, 

ikx k'^ Rii + kl R12 dy + ikx dyy, F^^ = 2 {ikx R12 dyy - kl Rn dy) . 
The operators appearing in the evolution equations ( |3.7[ ) are given by 

F2I + F2't;C^i + F2r]^r]l ^22 + F2't;C^2 + ^2r]^r]2 

^uv^vl + C^^Ct^i C^^C^2 + ^ur]^r]2 

^vl ^v2 
^wv^vl H~ ^wrj^rjl ^wv^v2 H~ ^wr)^r)2 

D = 



B 



Fi^D^ + Fi^D^ 

F2v'Dy + F277D^ 



{ikx / k'^) dy , ^uT] — ikz/k'^-) — {ikz/k'^)dy, C'ujyi = ikx/k'^. 



Appendix B. Explicit scaling of worst-case amplification of 
streamwise-constant velocity fiuctuations 

In this section, we discuss how to obtain exphcit expressions for the worst-case amph- 
fication from the wah-normal and spanwise forces to the streamwise velocity fluctuations 
in the limit of infinitely large We or infinitely large L. We note that derivation of the 
analytical expressions for Gu2 and Gus is more challenging because the worst-case am- 
plification of u arising from d2 and ds depends on both We and L. However, since our 
computations presented in § |5.1| demonstrate that the worst-case amplification from d2 
and ds to u takes place at low temporal frequencies, the essential features can be captured 
by analyzing the corresponding frequency responses at = 0; see figures [7(b)] and [7(^ 
For flows without temporal variations, we can obtain the following static-in-time expres- 
sions that relate the conformation tensor fluctuations with the streamwise velocity u and 
(?/,z)-plane streamfunction i/j (i.e., v = i/c^V^, w = —dy^l)) 

2We 2We We, ^ ^ 

r22 = -^ikzdy^l), rss = -^ikzdy^l), r2s = — YV^yy^^z)^^ (Bla) 

ris = {We/f) {U'r2s - {We/f) dyy^ + ikzu) , (B16) 
ri2 = {We/D (/7V22 - rii + {We/f) ikzdy^ + dyu) , (Blc) 

rn = {2WeyCo){{^We/f)ikzdy^ ^2dyu). {Bid) 
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The above expression (B 1 ) is obtained by taking the temporal Fourier transforms of (5.1a)- 
(5.16) and replacing the wall- normal velocity and vorticity with 

In the absence of streamwise forcing and streamwise variations, the static-in-time momen- 
tum equation (in the streamwise direction) provides a relation between the streamwise 
velocity and the streamwise components of the conformation tensor 



13 



13 f 



We 



dyri2 



We 



£2 



(B2) 



Substituting (B16) and (Blc) into (B2) yields a relation between streamwise velocity 
and the conformation tensor fluctuations in the wall-normal/spanwise plane, 



{l-p){dy{U'r22) + iA:.(t/V23)). 



(B3) 



Furthermore, we can obtain an expression relating u and ip by substituting (Bla 
into ( B 3 ) which yields 

Aix = - (We//) (1 -/3)ifc^ AV^. (B4) 
It can be shown that tjj is induced by the action of the wall-normal and spanwise forcing 



[ -ifc, dy ] 



d2 



(B5) 



Finally, the frequency response from d2 and da to at = is obtained by substitut- 
ing ( B 5 ) into ( B 4 ) which yields 

U^,{k,,0;l3,We,L) = -{We/f){l - /?)D„„ j = {2,3}, 

where 

D^2 = klA-^, D^3 = ik,A-^dy. 
The worst-case amplification of u caused by d2 and ds can be reliably approximated by 

Guj{h;p,We,L) ^ a^x (Hn, (fe., 0; /3, I^e, L)) 

= {We/f)' (l-/3)VLx(D..) 

= {N,/2){l-pfguj{k,), j = {2,3}, 

where the functions guj with j = {2,3} quantify the spanwise frequency responses from 
6^2 and ds to u. In the Oldroyd-B limit (i.e., as L ^ oo), the first normal stress difference 
Ni ^ 2 We^ and the function Guj is given by 

lim Guj{h;/3,We,L) = We^ {I - (if guj{k,). 

L^OO 

On the other hand, in the limit of infinitely large We, the first normal stress difference 
Ni and the function Guj is given by 



hm Guj{h;p,We,L) = {L^/2){1 



9uj {kz)' 
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Appendix C. The frequency response operators from body forces to 
polymer stresses in streamwise-constant Couette flow 
of FENE-CR fluids 



The streamwise-constant frequency response operators in (6.2) that map different forc- 



ing components to the polymer stress fluctuations are given by 

2/ 2/ 

Gij{k^,uj;/3,We,L) = . ^ , j dyT>yj, Gsj {k^, uj; We, L) = ^ 

lujp + / lujp + / 

G2j{k,,uj;p,We,L) = {f/{iujp^f)){i/k,){dyy^kl)B,j, 
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Wef 
(ic^/3 + /)(ic^ + /) 
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^ TT"; J ■E'uu D77I, 
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H^e(ic. + 3/) , E-iE D 
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